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As a model of a novel superconductor Na^CoC^-yHbO, a single-band t-J model on a triangular 
lattice is studied, using a variational Monte Carlo method. We calculate the energies of various 
superconducting (SC) states, changing the doping rate 5 and sign of t for small J/\t\. Symmetries 
of s, d, and d+id (p+ip and /) waves are taken up as candidates for singlet (triplet) pairing. In 
addition, the possibility of Nagaoka ferromagnetism and inhomogeneous phases is considered. It 
is revealed that, among the SC states, the d+id wave always has the lowest energy, which result 
supports previous mean-field studies. There is no possibility of triplet pairing, although the /-wave 
state becomes stable against a normal state in a special case (<5 = 0.5 and t < 0). For t < 0, the 
complete ferromagnetic state is dominant in a wide range of 5 and J/\t\, which covers the realistic 
parameter region of superconductivity. 



I. INTRODUCTION 

The recent discovery of superconductivity [l| in 
Nao.35CoO2d.3H2O has aroused an active interest of 
solid state physicists, since this compound is regarded 
as a strongly correlated superconductor on a triangular 
lattice. For such magnetically frustrated lattices, novel 
superconducting (SC) features are expected, as compared 
to the celebrated cuprates on bipartite lattices. In par- 
ticular, it is primarily important to determine the sym- 
metry of its pairing potential. So far, there have been a 
number of experimental efforts trying to fix the pairing 
symmetry, e.g. measurements of nuclear spin lattice re- 
laxation rate 2\ 0,0,0, Knight Shift 0,0 by NMR, and 
by fiSR At present, however, the results of these 

researches are not necessarily consistent to each another; 
the pairing symmetry has not yet been established. 

As for theory, although some studies regard multi- 
band effects as important, 0, 0, [^} many treat this 
issue with reduced single-band models on a triangu- 
lar lattice. The single-band studies are roughly clas- 
sified into two categories, namely ones applying weak- 
coupling techniques to Hubbard-type models, like pertur- 
bation expansions a random-phase approximation 
h2| and a fluctuation exchange approximation (FLEX) 
[l3|. and ones starting with t- J- ty pe models in the 
strong-coupling limit. [3, [H, [H, H3, HU Although the 
weak-coupling studies took mutually different approaches 
for different models, they arrived at a conclusion that 
/-wave spin-triplet symmetries are competitive with or 
even dominant over d-wave spin-singlet symmetries. On 
the other hand for i-J-type models, it is unanimously 
shown that d x 2 _ y 2 +id xy -w&ve symmetries are widely sta- 



ble against other symmetries as well as metallic states, 
using slave-boson mean- field approximations [T3. Il5l fl6 | . 
which barely respect the local constraint of no double 
occupation, and a Gutzwiller approximation^^, which 
takes account of the constraint on an average. Thus, the 
mean- field results of the strong-coupling theories disagree 
with the results of the weak-coupling ones, in contrast 
with the case of a square lattice, where the conclusion of 
the d x 2_ v 2 wave has been drawn in both theories. 



In this paper, we study the t-J model on a trian- 
gular lattice, using a variational Monte Carlo (VMC) 
method, which exactly treats the local constraint of 
no double occupation. Our main purposes are (i) to 
determine the pairing symmetry of superconductivity 
within the t-J model with this more accurate method, 
and (ii) to check the possibility of spin-triplet super- 
conductivity in this strong-coupling model. Besides 
the SC states, (iii) we consider the stability of Na- 
gaoka ferromagnetism [l9l I20L l2ll |22| and inhomogeneous 
phases. Based on these calculations, we argue the pro- 
priety of modeling for cobaltates. 



The organization of this paper is as follows: In §2, the 
formulation used is described. Sections 3 and 4 are as- 
signed to the results and discussions on the stability of 
the singlet and triplet SC states, respectively. In §5, we 
treat ferromagnetism and phase separation. In §6, we re- 
capitulate our results and address the relevance to cobal- 
tates. A part of the present results has been reported. [T^ 
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FIG. 1: (a) Upper panel: Schematic plot of the band struc- 
ture for Nao.5Co02, following an LDA calculation |24|. Lower 
panels: Dispersion of single-band tight binding models on a 
triangular lattice for t < and t > 0. (b) Hole Fermi surfaces 
of the CoO planes in Nao.sCoO^, due to the same LDA. 



II. 



FORMULATION 



eracy at e = ep is serious, we put t' — t and J' = J. 
Here, S indicates the hole concentration: 5 = |1 — n\ 
with n = N e /N s (N e : electron number, N s : site number) 
. This parameter setting makes no difference with the 
isotropic triangular lattice. 

Let us make a remark on the sign of t. In the lower 
panels of Fig.QJa), the tight-binding bands of the present 
model, 



Generally, the 3d orbitals are split into ti g and e g 
orbitals by an octahedral crystal field. In this ma- 
terial, the t2 g orbitals of Co are further split into 
the levels of a\ g [= (d xy + d yz + d Z3; )/v / 3] and e' g 



£k 



-2t (cos k x + cos k y ) — 2t' cos (k x + k y ) 



(2) 



3 d yz + e ±l s d 



yd X y 

oxygen vertical to the C0O2 layers. Due to an LDA cal- 
culation for Nao.5Co02[2J|, the bands near the Fermi 
surface [Fig. ^b)] are involved in a Co t 2g manifold 
[Fig.^a)]. Since Nao.35CoO2T.3H2O has the same struc- 
ture as Nao.5Co02, except for the intercalation of wa- 
ter molecules and the sodium concentration, one may 
assume the conduction bands for Nao.35CoO2T.3H2O 
are similarly formed from the t% g manifold. The band 
which forms the large Fermi surface around the T point 
(Fig.^b)) has dominant a\ g character, whereas the band 
which forms the small pocket Fermi surfaces on the T-K 
lines have mixed a\ g and e' g character. Thus, we call the 
band indicated by dashed (thick solid) lines in Fig. da) 
the aig (e' g ) band. Many previous studies have assumed 
that Nao.35CoO2T.3H2O is basically described by single- 
band models of the a\ g or e' g band, as mentioned. 

Along this line, we consider a single-band t-J model 
on an extended square lattice, 

U = H t + Hj = - tijVoicbcj* + H.c.)P G 

+ Y / J ij (s i S j --n i n j ), (1) 

where Vq = Y[ t (1 — n i] n ii)- The sum of site pairs 
is restricted as follows: — t and Jy = J(> 0) for the 
nearest-neighbor pairs, tjj = t' and = J'(> 0) for 
the next-nearest-neighbor pairs only in one diagonal di- 
rection, and tij = Jij = otherwise (see Fig. Except 
for the case of t < and i5 ~ 0.5 (§4), where the degen- 



with t' = t are depicted for t > and t < 0. Compared 
with the LDA band structure (upper panel of Fig.^a)), 
it seems appropriate that the band eq. @ with t < 



J /V3 [23, by a shift of [$ > 0] should be connected to the a lg (dashed lines) 



(bold lines)] band. In Na2;Co02-j/H20, the electron den- 
sity is larger than 1, i.e. electron-doped case. To treat 
electron-doped cases in the framework of t- J-type mod- 
els, namely in the electron density of < n < 1 without 
double occupation, we take advantage of a particle-hole 
transformation. Thereby, an electron-doped case with 
t < 0, which corresponds to the a\ g band (lower left panel 
of Fig.^a)), is mapped to a hole-doped case with t > in 
t-J model of eq. (JTJ. We use this sign convention through- 
out this paper. In fact, the t < band, which corresponds 
to the e' g band, has a van Hove singularity at S = 0.5, as 
the e' g band does near the K points. So far, the majority 
of theoretical studies[H El E E IH Q have con- 
sidered that the large Fermi surface around the T point 
mainly contributes, whereas some papers have regarded 
the isolated pocket Fermi surfaces or the van Hove sin- 
gularities as important for superconductivity, jjj E, Q 
In this paper, we consider both cases, t < and t > 0. 

Since the cobalt atoms in the present compound are 
connected through edge-shared oxygen cages, the effec- 
tive value of J/\t\ is expected to be fairly small. Accord- 
ingly, we principally discuss the region of small J/\t\ in 
this paper. 

We move on to the method used. To consider strongly 
correlated systems, more accurate treatment of the local 
constraint is necessary beyond the mean field level. For 
this purpose, a variational Monte Carlo (VMC) method 
ptl Ell l28j is appropriate. In this method, variational 
expectation values of many-body wave functions are nu- 
merically estimated with a Monte Carlo procedure, in 
which the local constraint, e.g. Vq in eq. (IJ, is exactly 
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TABLE I: Ak for various pairing symmetries. 



Symmetry 


A k 


s 


A 


ext. s (1) 


A (cos k x + cos k y ) 


ext. s (2) 


A[cos k x + cos(fc^ + k y ) + cos ky] 


d 


A(cos k x — cos k y ) 


d+id 


A cos k x + e l ~ cos(fe + k y ) + e % ~ cos k y 


p+ip 


A sin k x + e l2 $~ sm{k x + k v ) + e li S~ sin k y 


f 


A [sin k x — sin(k x + k y ) + sin k y ] 



satisfied. Actually, it has been shown that the VMC 
method is effective for a variety of analyses of the t-J- 
type models on square lattices without (or with small) 
frustration [2^, At least, we believe this method 

is reliable to discuss the mutual stability among various 
pairing symmetries. In this paper, we apply this method 
to the model eq. (JIJ to study the stability of various SC 
states. 

For i-J-type models, a Gutzwiller-type variational 
wave function is a good starting point, 13(1 HOJ un- 
like Hubbard- type models. [s3 For the normal state, we 
use the original Gutzwiller wave function ty n = Vg^f 
(<1>f: Fermi sea).[33| A SC state with a fixed electron 
number is given as, 

|* sc ) = P G |BCS) = P G f^ Vk 4 T c t _ k ^ |0), (3) 
with 

_ Mk _ A k 

k v k e k + V £ l + A k ' 

Anisotropy of the pairing symmetry is introduced in A k . 
We consider s, extended-s (two kinds), d and d+id waves 
as candidates for singlet pairing, and p+ip and / waves 
for triplet pairing. They are summarized in Table I. Sym- 
metries of s, p+ip and d+id waves have no line node in 
A k , whereas d and / waves have line nodes. The time- 
reversal symmetry is broken for d+id and p+ip waves. 
In the present calculation, A is the only variational pa- 
rameter, which is connected to the SC gap, in particular 
for a high-doping region, but is not the gap itself in con- 
trast with weak-coupling theories. Except for half filling, 
we may consider ty sc with finite A is a SC state, while at 
half filling \& sc is an insulating RVB state. [35| We fix // at 
Ho, the chemical potential for the noninteracting system, 
because fj, dependence of E to t is small, and \x does not 
seem an essential variable in a fixed- N e formulation. 

Since the number of variational parameters is at most 
one in this study, an usual VMC procedure |28l |34| suf- 
ficiently works. To obtain an accurate energy difference 
(condensation energy), we have taken a sufficient num- 
ber of samples to reduce the statistical errors, and kept 
the sampling interval long enough to ensure statistical in- 
dependence of the samples. Here, we collect samples as 



many as 10 5 -10 7 , which yield the error in energy of ap- 
proximately 10~ 4 \t\. We use the systems of N s = L x L 
(L=4-18) with the periodic-antiperiodic boundary con- 
ditions, and the electron densities satisfying the closed 
shell condition. 



III. STABILITY OF SINGLET PAIRINGS 

In this section, we discuss the symmetries of singlet 
pairing, s, extended-s, d and d+id waves. 

To begin with, we consider s-type pairings. Figure |3| 
shows the hopping energy E t = (/tit) and exchange en- 
ergy Ej = (Hj) for the homogeneous s-wave symmetry, 
as a function of A. For both t < and t > 0, each of 
E t /\t\ and Ej/J is a monotonically increasing function 
of A, indicating the energy minimum is at A = for 
any value of J/t; the s-wave state is not stable against 
the normal state. We also confirmed for various values 
of S that this situation does not change. This result is 
natural for the homogeneous s wave, which favors on- 
site pairing. As for the two kinds of extended-s waves, 
types (1) and (2), although E t /\t\ is a rapidly increasing 
function of A, Ej/J has a shallow minimum at a finite 
value of A. This means these symmetries become stable 
for finite values of J/t, unlike the homogeneous s wave. 
We confirmed, however, that the energy minimum is at 
A = for a small value of J/|t|.[3q] We thus conclude 
that there is no possibility of the s and extended-s waves 
realized in the regime of our concern. We will not refer 
to these states, henceforth. 

Next, we turn to the d- and d+id- wave pairings. Fig- 
ures 01 and El show E t /\t\ and Ej/J, respectively, of the 
d+id wave as a function of A/|i|. In this case, E t /\t\ is 
again a monotonically increasing function of A, whereas 
Ej/J conspicuously decreases and has a minimum at 
considerably large values of A. Thereby, the d+id state 
becomes stable even for considerably small value of J/\t\ 
(e.g. J/t ;> 0.05 for 6 = 0.04 and t > 0). Similar behav- 
iors of E t and Ej can be seen for the d-wave state, as 
seen in Figs. and [7\ From these figures, we find that 
the d+id and d waves have similar energies. To see actual 
energy reduction, we plot in Fig. OH the total energy i? to t 
(= E t + Ej) for J/\t\ = 0.3 and 6 = 0.04 and 0.12, as an 
example. For both d and d+id waves, -Etot has a clear 
minimum at large values of A for both signs of t; the SC 
states of d and d+id waves are appreciably stable near 
half filling. For S ~ 0, the minimum of the d+id wave is 
slightly lower than that of the d wave.[38j 

To consider the 5 dependence of energy, we define a 
condensation energy, E c , as the difference of .Etot be- 
tween the minimized superconducting state and the nor- 
mal state, 

E c =E(A min )-E(0). (4) 

For E c < 0, the superconducting state is stable against 
the normal state. Plotted in Fig. is E c for the d and 
d+id waves; \E C \ is a rapidly decreasing function of S. 
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FIG. 3: (a) Hopping energy and (b) exchange energy of the 
homogeneous s-wave state as a function of A/|f| for 8 = 0.04. 
The two cases, t > and t < 0, are simultaneously shown. We 
take the origin of each energy as the value for A == 0, namely 
the normal state, Et(0) and Ej(0); similar metrics are used 
in some figures following. The system used has 10 x 10 sites. 



Near half filling, the d+id-wave state is more stable than 
the d-w&ve state for both signs of t, probably because the 
symmetry of A& of the former is more favorable to the 
band symmetry, which is not far from isotropic. The dif- 
ference of E c , however, again becomes indistinguishable, 
as 8 further increases, and \E C \ itself becomes small. 

Note that E c of the d and d+id waves remains finite 
up to 8 = 0.24 for t > 0, and 8 = 0.60 for t < 0. Namely, 
the range of 8 where the d+id-w&ve state is dominant 
is much wider for t < than for t > 0. This behavior 
agrees with the Gutzwiller approximation even quan- 
titatively. This difference is understood as follows: For 
t < 0, the density of states p(e) increases as 8 increases, 
and has a van Hove singularity at 8 = 0.5, while, for 
t > 0, p(e) is a monotonically decreasing function of 8 
until the band bottom. This difference yields different 
behaviors particularly in Ej, as seen in Figs. and 
For t > 0, the minimum in Ej smoothly increases as 8 
increases for both d and d+id waves. On the other hand, 
for t < 0, the minimum in Ej once decreases as 8 in- 
crease, and has the lowest value at 8 ~ 0.36 for the d+id 
wave [Fig. Efb)]; for the d wave, the minimum in Ej is 
stationary for 0.04 <, 8 £ 0.20 [Fig.^b)] . For both waves, 
the decrease in Ej is still conspicuous for 8 as large as 



FIG. 4: Hopping energy of the d-Hd-wave symmetry for sev- 
eral values of the hole density 5. (a) t > and (b) t < 0. The 
system size is 10 x 10. Statistical errors are smaller than the 
symbol sizes. 

0.6. For this reason, the stability of the SC is enhanced 
particularly near the van Hove singularity point; we will 
return to this point in i jlVI 

The above VMC results that the d+id wave is domi- 
nant qualitatively support the conclusions derived from 
mean- field approximations 13 El an d a Gutzwiller 
approximation. |17| For comparison, we plot in Fig llOl 
.Etot(Amin) of the d+id-w&ve state by the present cal- 
culations and that of the Gutzwiller approximation. Al- 
though _B to t(A m i n ) by VMC is somewhat lower than that 
by Gutzwiller approximation, the two curves exhibit sim- 
ilar doping dependence. 



IV. POSSIBILITY OF TRIPLET PAIRINGS 

Ferromagnetic fluctuation, favored by triplet supercon- 
ductivity, has been predicted by an LDA calculation, [39| 
and actually observed in an NMR-Ti measurement for a 
superconducting sample. 4] In addition, many theories in 
the weak-coupling region have made conclusions of prob- 
able triplet superconductivity. In view of these sugges- 
tions, we study the stability of triplet symmetries, p+ip 
and / waves, in this section. 

For triplet symmetries, the same VMC procedure is 
carried out in the S z = sector. Having calculated ener- 
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FIG. 5: Exchange energy of the d+id-wave symmetry for sev- 
eral values of the hole density 8. (a) t > and (b) t < 0. The 
inset in (b) is the magnification of the minimum area. The 
system is the same as FigHI 



gies for various parameter values, it is found that Ej is 
always a monotonically increasing function of A. This is 
because the exchange term Hj/J with J > in cq. is 
obviously unfavorable to a triplet state. Moreover, E t /\t\ 
is also a monotonically increasing functions of A, except 
for the case of 8 ~ 0.5 and t < 0. Namely, the triplet 
superconducting states are always unstable against the 
normal state, excluding this special condition. In the 
following, therefore, we concentrate on this special case, 
t < and 8 = 0.5, at which the Fermi energy is situated 
just at the van Hove singularity point. 

At this point, severe level degeneracy arises due to the 
divergence of p(£f)- As a result, for finite systems, the 
electron densities where the closed-shell condition is sat- 
isfied disappear near 5 — 0.5. For example, this range 
is 0.4 < S < 0.6 for L = 10, as seen in FigEJb). To 
avoid this difficulty, we slightly change the value of if /t 
from 1 to e.g. 0.999999, which lifts the degeneracy and 
make it possible to study 8 = 0.5 without breaking the 
closed-shell condition. |4£| 

In Fig. II If . Ej is shown for four symmetries. As 
mentioned above, Jij works disadvantageously to the 
triplet pairing symmetries; Ej for the p+ip and / waves 
is again a monotonically increasing function of A, in con- 
trast with the behavior of the singlet pairing symmetries. 
Thus, the triplet SC states are never stabilized when J/\t\ 



FIG. 6: Hopping energy of the d-wave symmetry for several 
values of the hole density 8. (a) t > and (b) t < 0. The 
system size is 10 x 10. 

increases. In Fig. lllf aV Et is depicted for the same sym- 
metries. Note that Et for the / and p+ip waves (solid 
symbols) has a shallow minimum at a finite value of A. 
From this result, we find that the /-wave state becomes 
stable against the normal state (A = 0) for small values 
o£J/\t\ (J/\t\ ^0.18). Compared with the singlet pairing 
symmetries, however, the decrease in E t for the triplet 
pairings is much smaller than that for the d-\-id and d 
waves. Moreover, the singlet SC states gain energies in 
Ej for finite values of J/t. We add the data of the d 
and d+id waves at this density (J/\t\ — 0.3) to Fig.|§fb); 
we find that the van Hove singularity appreciably con- 
tributes to superconductivity. Nonetheless, we will see 
that even this d wave is defeated by ferromagnetism for 
small J/\t\. 

To summarize, triplet superconductivity is not realized 
in the triangular t-J model. 



V. FERROMAGNETISM AND PHASE 
SEPARATION 

In this section, we discuss ferromagnetic states and 
inhomogeneous phases, which are inseparable from the 
t-J model. 

Let us begin with the Nagaoka ferromagnetism. [l9| 
Since Nagaoka's theorem (in the limit of J/\t\, 8 — > 0) 
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FIG. 7: Exchange energy of the d-wave symmetry for several 
values of the hole density S. (a) t > and (b) t < 0. The 
inset in (b) is the magnification of the minimum area. The 
system is the same as Figltll 



FIG. 8: Total energy of the d- (solid symbols) and d+id-wave 
(open symbols) states for 5 = 0.04 and 0.12. The system is of 
10 x 10, and J/\t\ = 0.3. (a) t > and (b) t < 0. 



holds for t < in a triangular lattice, a finite range of 
ferromagnetism is expected near this limit. Actually, it 
was shown using a variation scheme for the Hubbard 
model on a triangular lattice with t < that a Na- 
gaoka itinerant ferromagnetism is widely stable. [2(j For 
the t-J model with t < 0, recent studies using a high- 
temperature expansion[2lj showed the Nagaoka region 
to be < S ^ 0.85 and J/t ^ 0.5, which is quite wider 
than that for a square lattice (possibly a partial ferro- 
magnetism), namely, < 6 £ 0.4 and J/t £ O.l.pll^ 
With these in mind, we look at our VMC results in the 
following. In Fig. 1121 we show E t for the complete ferro- 
magnetism (Nagaoka state), E { , and the d-Md-wave SC 
state, E? +ld . For t > 0, E { is always higher than E? +id ; 
thus the complete ferromagnetism never takes place, be- 
cause (TLj) = for the Nagaoka state, while E d+ld fur- 
ther decreases due to the contribution of Hj. Inversely, 
for t < 0, E f is appreciably smaller than E d+ld for most 
values of 5, thereby there exist a region for small J/\t\ 
where the Nagaoka state is more stable than the d+id- 
wave state. 

Comparing E l and E d+ld for finite values of J/t, we 
determined the phase boundary, as shown in Fig. I13f b). 
The range of ferromagnetism is very broad, namely 5 < 
0.96 and J/\t\ 0.7, which is quantitatively consistent 
with a recent result of high-temperature expansion. |2lj 



and is markedly wider than that for the square lattice 
estimated through the same VMC method. HllH Con- 
sequently, for t < 0, it is unlikely that superconductivity 
is realized in the region of small J/\t\ and intermediate S, 
even if we consider that the ferromagnetic region some- 
what shrinks by improving the trial function of supercon- 
ducting states. 

Now, we turn to the phase separation in the triangular 
lattice. Since the exchange term in eq. (JQ works as a 
kind of attractive potential, the t-J model necessarily 
brings about a phase separation for large values of J/\t\. 
For the one-dimensional lattice, the boundary between 
homogeneous and inhomogeneous phases extends in the 
direction of large J/t as 6 decreases ^lj, whereas for 
the square lattice, the boundary approaches the origin 

[(j/m) = (o,o)] mM 

Here, we estimate the boundary in two distinct ways, 
(i) We suppose that the inhomogeneous state phase sep- 
arates into the two hole densities, S = (half filling) and 
5 = 1 (empty) , and the energy of this state is estimated 
as E PS = (1 — S)E AF , where we adopt, as E AF , a pre- 
cise value obtained by a quantum Monte Carlo method 
for the Heisenberg model, (43 in which the antiferromag- 
netic (AF) state with the 120° structure is realized. The 
boundary is determined by comparing E PS and E d+ld . 
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FIG. 9: Condensation energy of the d- and d+id-wave states 
as a function of doping rate S for J/\t\ = 0.3. (a) t < and (b) 
i > 0. In (b), the data for 8 — 0.5 (the van Hove singularity 
point) is obtained under a different condition; for details, see 
ill VI Note that the range of abscissa is twice larger for t < 0. 
The Data for L — 10, 12 and 18 are simultaneously plotted; 
system-size dependence is not severe. 



Note that the reliability of this estimation increases as 
6 increases, but the boundary necessarily converges to 
J/\t\ — as 8 — > 0. 0| (ii) Instead, we estimate the 
boundary near half filling, using the condition of intrin- 
sic stability, d 2 E/d5 2 > 0, for the d+id-wave state. Since 
we actually substitute finite differences for d 2 E/dS 2 , the 
critical value of J/\t\ at half filling for L — > oo is likely to 
be smaller. 

The results of (i) and (ii) are summarized in Fig. ED 
For t < 0, the critical value J c /\t\ is relatively small and 
barely depends on 8 for large 5. For t > 0, J c /\t\ is 
fairly large for large S, and gradually decreases as S de- 
creases. For both signs of t, J c /\t\ abruptly decreases as 
S approaches half filling; this tendency is analogous to 
the square lattice case. Although the possibility of phase 
separation still remains for a realistic value of J/\t\ near 
half filling, this is not the case with intermediate and 
large values of 8. 
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FIG. 10: Comparison of £ , tot(A m i 11 ) of the d+id-wave pairin, 
between the present VMC and a Gutzwiller approximation^ 
for t < and J/\t\ = 0.3. 



VI. 



SUMMARY AND DISCUSSIONS 



In this paper, we have studied the single-band t-J 
model on a triangular lattice, based on variational Monte 
Carlo (VMC) calculations, with Nao.35CoO2T.3H2O in 
mind. The present study, which treats the local correla- 
tion precisely, is more reliable than the previous mean- 
field-type calculations. We have compared the total en- 
ergy for various pairing symmetries in a small- J/\t\ re- 
gion, changing the doping rate 8 and the sign of t. We 
consider s, extended-s, d, and d+id waves for singlet pair- 
ing, and p+ip and / waves for triplet pairing. Further- 
more, we have studied the stability of Nagaoka ferro- 
magnetism and inhomogeneous states. Main results are 
summarized as follows. 

(1) For t > 0, the d+id-wave state, which breaks time 
reversal symmetry, is dominant for < 8 ^ 0.24 (J/t = 
0.3), but the simple d-wave state has a quite close energy 
to that of the d+id wave for a similar range of 8. The 
stability of the superconducting state is rapidly lost, as 
the carrier density goes away from half filling. For 8 ~ 0, 
phase separation may take place. See Fig.^fa). 

(2) For t < 0, in a wide parameter range of our con- 
cern (0 < 8 £ 0.96, < J/t <, 0.7), the Nagaoka ferro- 
magnetic state defeats the d+id-wave state, which is the 
most stable among the superconducting states studied. 
See Fig. Eft)). 

(3) Triplet pairing is not realized in the t-J model, 
although an /-wave state becomes slightly more stable 
than the normal state in the special case of 8 = 0.5 and 
t < 0, where the Fermi surface overlaps the van Hove 
singularity points. 

These results broadly support the previous mean- 
field- type studies. [bH ITa. Hil Il7| and is consistent with 
the high-temperature expansion studies. Now, let us 
check these results for the triangular t-J model in the 
light of a model of Nao.35CoO2-l.3H2O. Within the t-J 
model with small J/\t\, the range of steady supercon- 
ductivity is limited to < 8 ^ 0.24 with t > 0, and 
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AM 



FIG. 11: Comparison of (a) hopping energy and (b) ex- 
change energy among the p+ip, f, d+id and d waves for 
t'/t = 0.999999 (t < 0) and 8 = 0.5 (L = 12). Here, we 
assume J' /J = (t'/t) 2 . The inset in (b) shows the behavior 
of the d and d+id waves in a wider range. 



the pairing symmetry is the d+id (or possibly pure d) 
wave. In experiment, although the pairing symmetry of 
Nao.35CoO2T.3H2O has notyet established, it is likely 
that its gap has a line nodeJ&Q and does not break the 
time-reversal symmetry. 0, Q Furthermore, the carrier 
density at which superconductivity appears is limited to 
a narrow range near S = 0.3. Therefore, it is difficult 
to conclude that the t-J model successfully describes the 
properties of Nao.35CoO2T.3H2O in the present exper- 
imental situation. Provided that this discordance orig- 
inates in the theoretical side, we discuss some related 
issues briefly in the remainder. 

As mentioned in §1, the results of triplet, especially /- 
wave, pairing were obtained by applying weak-coupling 
approaches to Hubbard- type models. [ill Il2l| Among 

them, Kuroki et al. paid attention to the role of e' g 

band|13|. which broadly corresponds to the t < band. 
They showed using a FLEX approximation that the /- 
wave state is the most stable. Thus, it is interesting to 
confirm whether the stabilization of the / wave found for 
S ~ 0.5 and t < in the t-J model is connected to this 
weak-coupling / wave in the Hubbard model. Actually, 
the Hubbard-type model became tractable by a VMC 
method due to a recent improvement of the variational 



FIG. 12: Comparison of the transfer energy between the op- 
timized VMC values (J/t = 0) of the d+id-wave SC state and 
those of the complete ferromagnetic state (E 1 ), estimated an- 
alytically through a spinless Fermion model, (a) t > 0, and 
(b) t < 0. 



wave functions. |4fij 

In this paper, we have used a single-band model as a 
first step to study the Co-based oxide. However, there 
exists a possibility that the multiband effect plays an 
essential role for superconductivity. H, 0] Actually, re- 
cent FLEX and perturbation calculations performed 
on multi-band Hubbard models yielded results of dom- 
inant / and p waves. Here, the e' g band again makes 
a key contribution to superconductivity. To study this 
effect for large U jt is also an interesting remaining issue. 

In the connection of triplet pairing, we finally remark 
that the pairing symmetry of Nao.35CoO2T.3H2O can be 
determined by the phase sensitive experiments like tun- 
neling effect and Jose phs on effect through the Andreev 
bound states. H3, EEHS E3, H3 In particular, to iden- 
tify the triplet pairing, it is effective to use the anomalous 
proximity effect predicted very recently. p5^ 
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FIG. 13: Phase diagram in the S-J/\t\ space constructed from 
the present calculations, for (a) t > and (b) t < 0. The 
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phase separation estimated by the way (ii). Note that the 
system-size dependence is insignificant. 
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